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INTRODUCTION 
A  Definition  of  Dynamic  Stability 

In  the  beginning  of  time  on  this  earth  unstability  in  motion  was 
observable  when  the  first  leaf  fell  from  the  first  plant.  As  the  leaf 
descended  teward  the  earth  it  either  fluttered,  oscillated,  revolved, 
spiraled,  or  combined  all  or  part  of  these  motions  in  its  erratic  descent. 
This  erratic  behavior  was  indicative  of  the  non-uniform  velocity  of  the 
leaf  along  a  vertical  to  the  earth.  Therefore  it  was  exhibiting  unstable 
motion,  or  dynamic  unstability. 

If  the  first  man  had  been  able  to  measure  the  height  of  the  leaf  at 
equal  increments  of  time,  dynamic  stability  might  have  been  acknowledged 
millions  of  years  ago.  Initially  the  first  man  night  have  noticed  that 
the  measured  heights  of  the  leaf  did  not  decrease  the  same  amount  over 
each  equal  increment  of  time.  This  indicated  non-uniform  velocity  in 
the  vertical  direction.  The  first  man  might  then  have  compared  these 
measured  heights  with  the  comparable  heights  of  a  leaf  imagined  as 
descending  toward  the  earth  with  uniform  velocity  in  the  vertical 
direction.  From  comparing  each  actual  height  with  each  imagined  height 
for  the  same  time  of  descent  he  would  have  found  a  difference.  By  then 
calculating  the  mean  of  the  summation  of  these  absolute  differences  he 
might  have  called  the  resultant  value  dynamic  unstability.  However,  he 
didn't.  Dynamic  unstability  had  always  existed,  but  it  was  necessary 
for  man  to  measure  It  to  verify  its  existence. 

Simply,  dynamic  unstability  was  the  mean  deviation  of  the  actual 
path  of  the  leaf  from  some  imaginary  path,  described  in  this  study  by 
a  least  squares  line  through  the  actual  path.  For  this  imaginary  path 


no  assumptions  were  made  as  to  the  uniformity  of  the  leaf's  velocity 
vectors  horizontal  to  the  earth's  surface.  It  was  possible  for  the 
leaf  to  follow  a  curvelinear  path  in  its  actual  descent  and  exhibit  no 
dynamic  unstability  as  long  as  it  descended  equal  increments  of  vertical 
distance  in  equal  increments  of  universal  time. 

The  difference,  or  deviation,  of  the  actual  path  from  an  imaginary 
path  had  units  of  length.  Therefore,  dynamic  unstability  had  units  of 
length.  As  the  phraseology  might  indicate,  dynamic  stability  was  defined 
as  the  reciprocal,  or  opposite,  of  dynamic  unstability.  Therefore,  its 
units  were  reciprocal  length.  In  order  to  adequately  describe  the  meaning 
of  dvnamic  stability  it  was  necessary  to  define  it  in  mathematical  terms; 
this  was  presented  in  the  Mathematical  Definition  of  Dynamic  Stability. 

Many  investigators,  including  tfadell  (8),  Pettyjohn  (6),  Squires  (7), 
Wetherall  (9)5  and  Becker  (l),  had  observed  this  behavior  in  most  des- 
cribable  sizes  and  shapes  of  particles,  but  carried  the  investigation 
little  further  than  recording  these  observations  and,  in  some  cases, 
comparing  them  with  conventional  concepts  of  particle  behavior  in  fluids. 
No  attempt  had  been  made  to  numerically  evaluate  or  to  define  the  un- 
stability of  movement  of  particles  in  fluids. 

Disks  which  exhibit  infinite  stability  are  illustrated  in  PLATE 
I.  They  were  defined  as  a  stationary  disk,  floating  on  the  surface  of 
a  fluid,  with  no  velocity  vector  and  as  a  moving  disk,  descending  or 
ascending  in  the  fluid,  with  velocity  but  no  acceleration.  Conversely, 
infinite  unstability  was  postulated  as  the  motion  of  a  disk  which 
exhibited  constant  acceleration  through  infinite  time.  A  disk  accelerating 
and  decelerating  through  finite  time  had  a  finite  value  of  dynamic  stability 


and  unstability.  It  was  these  latter  descents  of  free-falling  disks  with 
which  this  study  was  concerned. 

Even  as  the  best  measurement  is  no  better  than  an  approximation  of 
the  true  value,  so  no  motion  is  absolutely  dynamically  stable.  Due  to 
limitations  on  man's  ability  to  measure  minute  differences,  it  is  possible 
to  measure  and  subsequently  calculate  a  descent  as  dynamically  stable, 
although  theoretically  a  dynamically  stable  descent  is  unobtainable.  As 
man's  abilities  to  measure  improve  that  which  was  previously  classified 
as  dynamically  stable  may  no  longer  be  so  classified. 

Objectives  of  Study 

The  original  objectives  of  this  study  weres 

1)  To  continue  work  initiated  by  Professor  R.  C.  Hall*  and  to 
refine  the  mathematical  definition  of  the  concept  of  free- 
falling  particle  behavior  in  an  infinite  expanse  of  fluid. 

2)  To  calculate  numerical  values  of  stability  from  actual 
descents  of  disks  of  varying  properties. 

Later  objectives  weres 

3)  To  observe  and  define  the  types  of  motion  of  free-falling 
disks  in  a  fluid  media. 

4)  To  correlate  these  values  of  dynamic  stability  with  the  types 
of  motion  observed. 

Objectives  3)   and  u)   were  planned  after  an  initial  study  of  the  types 
of  motions  of  disks  had  shown  that  some  of  the  disks  consistently  fell  in 
the  same  manner. 

♦Private  communications 


Review  of  Literature 

A  literature  search  revealed  that  dynamic  stability  had  aroused  the 
curiosity  of  many  investigators,,  but  that  this  effort  was  apparently  the 
first  actually  measuring  it.  The  literature  search  included  the  follow- 
ing materials 

1.  Chemical  Abstracts  dating  back  to  1930. 

2.  Physics  Abstracts  back  to  1950. 

3.  All  theses  and  abstracts  of  theses  located  in  the  Kansas  State 
University  Library. 

U.     Numerous  more  recent  publications  not  yet  abstracted. 

No  experimentation  resembling  dynamic  stability  studies  was  un- 
covered. However,  several  articles  discussing  the  appearance  of  motion 
in  particles  descending  through  a  fluid  media  were  available. 

From  an  article  by  Wadell  (8),  the  following  was  quoted; 

Franz  Schulze  stated  that  the  translation  velocity  of  a 
flat  stone  in  air  or  water  iepends  upon  its  position  in  respect 
to  the  direction  of  translation.  Jordan  noted  that  mineral 
particles  readily  took  a  position  with  their  broad  side  at  right 
angles  to  the  direction  of  fall.  Schmiedel  found  that  thin 
disks  at  very  low  Re  showed  no  tendency  to  take  any  particular 
settling  position  while  at  higher  Re  the  disks  tended  toward  a 
horizontal  position  regardless  of  the  position  at  the  start.  At 
still  higher  Re's  periodic  oscillations  set  in  and  the  settling 
took  place  in  a  zig-zag  path. 

On  the  study  of  spheres,  octahedrons,  cubes,  and  tetrahedrons 

Pettyjohn  and  Christiansen  (6),  stated* 

Wheras  in  the  streamline  flow  range  (Re<0. 05)  the 
particles  did  not  favor  any  particular  orientation  with  re- 
spect to  the  direction  of  motion,  at  Re'^lO,  the  tetrahedrons 
and  cubes  assumed  an  orientation  with  a  face  in  the  horizontal 
plane  or  perpendicular  to  the  direction  of  motion  and  the 
other  particles  showed  a  tendency  to  do  likewise  by  Re<-^20. 
This  orientation  was  maintained  up  to  Re  ==70  to  300  where  the 
particles  beginning  with  the  tetrahedrons,  'teetered'  or 


wobbled  and  eventually  spun  or  rolled  on  a  horizontal  axis 
and  followed  a  spiral  path  rather  than  a  straight  vertical 
path  in  their  descent. 

From  an  article  on  the  sedimentation  of  thin  disks  by  Squires  and 

Squires  (7),  it  was  stated  that  at  very  low  Reynolds  numbers,  in  the 

viscous  region  a  body,  such  as  a  disk,  possessing  three  perpendicular 

planes  has  no  tendency  to  assume  any  particular  orientation  during 

fall.  It  was  further  stated: 

A  disc  placed  with  its  flat  face  parallel  to  the  direction 
of  motion  of  the  fluid  will  maintain  this  orientation.  At  higher 
Reynolds  numbers,  in  the  turbulent  region,  the  eddies  set  up  by 
the  motion  of  the  disc  through  the  fluid  act  to  maintain  the 
particle  at  right  angles  to  the  direction  of  motion.  Therefore, 
a  disc  will  always  settle  in  a  horizontal  position  in  the  turbu- 
lent region,  and  its  position  when  in  the  viscous  region  will 
depend  upon  its  initial  inclination. 

Miyagi  (5),  discussed  the  motion  of  an  air  bubble  rising  in  water. 
The  course  of  the  bubble  was  three-dimensional;  it  was  observed  to  be 
a  helix  around  a  vertical.  The  major  axis  of  the  bubble  was  always 
perpendicular  to  its  course.  PLATE  XVIII,  Fig.  1,  based  on  Fig.  12  of 
Miyagi' s  article  illustrated  the  rise  of  an  air  bubble  in  only  two  of 
the  three  dimensions. 

Wetherall  (9),  was  able  to  photograph  the  track  of  a  flat  solid 
(a  small  crystal  of  A-110*)  as  it  descended  in  acidulated  water.  PLATE 
XVIII,  Fig.  2,  was  based  on  this  photograph. 

Becker  (l),  provided  the  following  information  in  a  table  on  the 
inertial  drag  characteristics  of  freely  oriented  bluff  bodies.  At 
Re' s  of  0.1  -  5.5  it  was  stated  that  all  orientations  were  stable  when 
there  were  three  or  more  perpendicular  axes  of  symmetry  in  the  descending 
disk.  At  Re's  of  5. 5  -  200  the  motion  was  stable  in  the  position  of 
maximum  drag.  For  Re's  of  200  -  500  the  motion  was  unpredictable. 


Disks  tended  to  wobble,  while  fuller  bluff  bodies  tended  to  rotate. 

Considering  the  results  of  this  literature  search  it  was  apparent 
that  this  study  initiated  by  Professor  R.  C.  Hall  of  Kansas  State 
University  in  1955  was  the  first  recognition  of  dynamic  stability  as 
a  property  of  particles  in  an  infinite  expanse  of  fluid. 

Illustration  of  Definition  of  Dynamic  Stability 


The  following  diagram,  PLATE  I,  illustrates  the  two  possible 
manners  in  which  a  disk  could  theoretically  exhibit  complete  dynamic 
stability. 


EXPLANATION  OF  PLATE  I 


Illustration  of  infinitely  stable  motion  of  disks, 
where  y  —  dynamic  unstability  =  0, 

^  =  dynamic  stability  =  CO   ,   and 

U  =  Telocity. 
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PROCEDURE  OF  MEASUREMENTS 


Outline  of  Procedure 


For  the  study  a  series  of  46  disks,  approximately  one  inch  in 
diameter  and  a  quarter  inch  thick,  were  constructed.   Disk  1  was 
lucite.   Disks  2  through  26  had  cores  of  stainless  steel  and  rims 
of  lucite.  Disk  27  was  all  stainless  steel.  Disks  28  through  46 
had  lucite  cores  and  stainless  steel  rims.  By  varying  the  size  of 
the  cores  the  disks  were  each  given  different  weights  and  moments 
of  inertia.  Only  part  of  the  disks  were  used,  due  to  the  deterior- 
ation of  the  rest. 

The  procedure  was  to  drop  each  disk  in  turn  through  a  liquid, 
take  motion  pictures  of  the  descent,  and  measure  on  a  projection  of 
each  frame  of  the  movie  the  disk's  relative  height  above  a  bottom 
reference  line  (PLATE  XII). 

Detailed  Procedure 


The  procedure  utilized  in  taking  the  movies  of  the  disk's  descent 
in  a  liquid  media  required  the  services  of  assistants  A  and  B  and 
photographer  P  (PLATES  II,  III,  and  IV).  For  the  descent  of  each 
disk  the  procedure  was  as  follows! 

1. )  As  assistant  B  changed  the  cards  designating  the  disk 
number,  the  drop  number,  and  the  liquid  temperature, 
assistant  A  placed  the  designated  disk  on  the  particle 
vacuum  release  mechanism  performing  operations  1,  2, 
and  3  (PLATE  V). 
2. )  B  then  turned  on  the  foreground  lights  illuminating  the 
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information  cards.  This  signaled  P  to  photograph  this 
information. 
3).  Having  built  up,  during  a  twenty  second  count,  sufficient 
vacuum  to  hold  the  disk  the  line  from  the  vacuum  pump  was 
closed  and  "bled1  by  A  in  performing  operations  U  and  5 
(PLATE  V)„  The  release  mechanism  holding  the  disk  was 
then  lowered  into  the  liquid „ 
4. )  At  the  end  of  fifty-five  seconds  of  a  second  count,  A 
notified  B  to  signal  P  by  means  of  a  buzzer  to  commence 
photographing,,  At  sixty  seconds  the  remaining  vacuum 
was  released  by  performance  of  operation  6  (PLATE  V), 
thus  allowing  the  disk  to  start  its  descent  through 
the  liquid. 
5. )  As  the  particle  reached  the  bottom  of  the  column,  B 

signaled  the  photographer  to  cease  operation. 
6. )  This  procedure  was  repeated  for  successive  disks. 
An  approximately  identical  method  of  release  for  the  disks  was  ob- 
tained by  means  of  the  particle  vacuum  release  mechanism.  By  allowing 
a  time  lapse,  the  twenty  and  sixty  second  counts,,  and  by  floating  Agile 
plastic  stars  (See  Table  of  Definitions,  Appendix)  on  the  fluid  surface, 
the  internal  fluid  agitation  was  minimized  to  a  satisfactory  degree. 

Following  the  development  of  the  film,  measurements  on  the  descent 
of  the  disk  were  started.  For  this  operation  assistants  A  and  B  were 
again  utilized.  The  assistants  were  selected  on  the  basis  of  the 
similarity  between  their  measurements  over  several  trials.  This  like- 
ness in  observation  of  the  measured  quantity  allowed  the  assistants  to 
be  rotated  between  jobs  within  any  set  of  measurements.  To  further 
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minimize  fatigue  and  monotony  a  third  assistant  C  was  used  when  avail- 
able. 

The  film  was  threaded  through  a  special  adaptor  (PLATE  XI)  on  the 
projector.  The  use  of  this  adaptor  largely  prevented  the  film  from 
shifting  or  buckling;,  as  well  as  adapting  the  16mm„  film  to  a  35mm. 
projector.  Prior  to  starting  the  measurements  the  projector  was 
allowed  to  run  twenty  to  thirty  minutes ,  thus  allowing  distortion 
caused  by  heat  from  the  projection  lamp  to  reach  near  equilibrium  in 
the  film.  The  focus  on  the  projector  was  earlier  set  at  an  optimum 
position;  this  was  determined  by  repeated  adjustments  of  the  focus 
until  both  assistants  were  visually  satisfied,, 

The  actual  measuring  operation  was  in  general  performed  as 
follows : 

1. )  After  the  projector  was  allowed  to  heat  for  twenty  to 
thirty  minutes  before  the  day's  run,  the  disk  number, 
liquid  temperature,  and  drop  number  were  identified  by 
assistant  B  from  an  identification  frame  (PLATE  XII, 
Fig.  l)  and  called  to  assistant  A  to  be  recorded. 
2. )  The  film  was  then  pulled  through  the  projector  by 
assistant  A  until  the  frame  was  found  showing  the 
disk  to  be  located  just  below  the  top  reference 
line.  This  frame  was  the  first  to  be  measured  for 
each  disk  descent. 
3.)  At  this  point  necessary  minor  adjustments  were  made 
in  the  frame  location  on  the  screen.  By  means  of  the 
turnbuckles  and  the  focus  on  the  projector  (PLATE  VI 
and  VII)  the  projected  picture  of  the  frame  was  re- 
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strieted  between  vertical  boundary  lines  drawn  on  the 
screen  (PLATE  VIV).  By  pulling  the  film  up  or  down  the 
image  of  the  top  reference  line  was  placed  in  approximate 
line  with  a  reference  mark  drawn  on  the  screen  (PLATE  VIV)„ 
These  adjustments  minimized  any  possible  error  in  measure- 
ments. These  errors  occurred  from  minute  differences 
between  each  frame's  position  in  the  projector  due  to 
the  slight  bending  of  the  film  and  the  minor  deviations 
of  the  film  from  the  vertical  in  passing  each  frame 
through  the  adapter  (PLATE  XI),, 

4„ )  After  these  preliminaries  were  finished  the  readings  on  a 
frame  were  started.  By  means  of  a  target  (PLATES  VI,  IX, 
and  X),  B  located  the  bottom  reference  line,  the  particle, 
and  the  top  reference  line  positions,  respectively.  These 
locations  were  in  turn  transmitted  by  a  system  of  pulleys 
to  a  stationary  measuring  device  (PLATE  VI  and  VIII)  where 
the  values  Lg„  Lp,  and  L^  were  read  and  recorded  by  A„ 
For  reasons  of  uniformity  in  the  tension  on  the  pulley 
wire  the  target  was  always  adjusted  downward  onto  the 
image  being  measured.  For  statistical  reasons  the  readings 
were  repeated  six  times  for  each  frame.  For  each  reading, 
by  means  of  a  signal  light,  B  notified  A  when  the  target 
was  set  and  by  means  of  a  buzzer,  A  notified  B  when  he 
had  completed  the  reading  (PLATE  VIII). 

5. )  After  completing  this  first  frame  a  second  frame  was  placed 
in  position.  For  each  frame,  after  the  initial  frame  of  a 
day's  measurements,  a  ten  to  fifteen  second  time  lapse  was 
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allowed  so  that  the  frame  had  sufficient  time  to  reach 
equilibrium  with  the  distortions  caused  by  the  heat  of 
the  projection  lamp  striking  the  film. 

6. )  The  same  procedure  was  employed  for  each  frame.  The 
measurements  for  any  one  disk  were  terminated  at  the 
frame  showing  the  particle  in  its  final  position  above 
the  bottom  reference  line. 

By  utilization  of  a  true  length  correction  factor  for  each  frame,  . 

100  .  changes  in  the  film  due  to  moderate  temperature  and  humidity 
Lp  _  Ll 

changes  were  prevented  from  causing  measuring  errors.  More  rapid  dis- 
tortion due  to  extreme  temperature  changes  in  the  film  as  caused  by  the 
heat  from  the  projection  lamp  were  minimized  by  short  time  lapses 
between  each  frame  read  in  which  equilibrium  was  reached  in  the  film. 
Fatigue  in  the  measuring  personnel  was  minimized  by  changes  in  duties 
every  half  hour  to  one  hour,  depending  on  the  weather  and  by  taking 
frequent  rest  breaks. 

Errors  in  measurement  caused  by  shifts  in  the  picture  projection 
due  to  extraneous  disturbances  were  avoided  by  the  replacement  of  such 
data  after  the  disturbance  had  ceased.  Minor  disturbances  and  vibrations 
were  alleviated  by  the  location  of  rubber  padding  under  the  projector 
platform  and  by  the  use  of  several  restraining  wires  around  the  projector 
(PLATE  VI  and  VIII). 

Diagrams  of  Equipment 

The  following  pages  present  diagrams  of  the  experimental  equipment 
and  illustrate  procedures. 


EXPLANATION  OF  PLATE  II 

Front  view  of  particle  drop  column  showing 
lighting  arrangement  and  personnel  placement 
for  photography  of  disk  descent  in  liquid 
column. 
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PLATE  II 


AGILE 

PLASTIC 

STARS 


eACKTROuND    LIGHTS 


EXPLANATION  OF  PLATE  III 

Top  view  of  particle  drop  column  showing  personnel, 
lighting  and  camera  positions  for  photography  of 
disk  descent  in  fluid  media. 
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PLATE  III 
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EXPLANATION  OF  PLATE  IV 

Schematic  of  arrangement  for  photography  of 
disk  descent  in  fluid  media  in  glass  column. 
/ 1     1 1  heavily  starched  white  sheet  as  back- 
ground of  drop  tank. 
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EXPLANATION  OF  PLATE  VII 

Front  view  of  projector  showing  location  of 
turnbuekles  for  right  and  left  adjustments 
of  film  projection  and  focus  for  refinement 
and  right  and  left  cant  of  projection. 
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PLATE  VII 
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EXPLANATION  OF  PLATE  VIII 

Front  view  of  measuring  device  showing 
assistant  A,  projectionist  and  reader- 
recorder  of  L^,  Lp,  and  Lgj  values. 
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PLATE  VIII 
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EXPLANATION  OF  PLATE  X 


Front  view  of  target,  showing  actual  application 
of  target  to  three  measured  elements  of  each  frame. 
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PLATE  X 
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EXPLANATION  OF  PLATE  XII 

Frames  from  movie  of  the  descent  of  a  disk. 
Fig.  1.  Identification  frame. 

///.  temperature  of  fluid  media 

>XS^\disk  number 

==.  drop  number 
Fig.  2.  Measuring  frame. 
Lj,  =   top  reference  line 
Lp  =   disk 
Lt  =   bottom  reference  line 
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MATHEMATICAL  DEFINITION  OP  DYNAMIC  STABILITY 
Initial  Calculations  on  Data  from  Representative  Frame  of  a  Descent 


6  =  frame  number. 

Lg  =  measured  height  of  top  reference  line,  relative  to  base  line. 

L_  •  measured  height  of  disk,  relative  to  base  line. 


L 
T.L. 


measured  height  of  bottom  reference  line,  relative  to  base 
line. 

-   100  =  factor  for  adjusting  relative  values  to  true  values. 
LH  "LL 


■ 

h.  =  (L„  -  L. )  (    100   )  =  true  height  of  disk  above  bottom 
L  -  LT    reference  line. 

(See  PLATE  XIII,  Fig.  1  for  physical  descriptions  of  L  ,  Lp,  and  LL) . 

For  first  reading  on  first  frame: 

&      Reading        L„             Lp             LL 

1        1        105.23         104.52         25.52 

LH  -  LL  -  79.71 

Lp  -  LL  =  79-00 

T.  L.  =    100    =    „,, 

LH  "  LL 

h.  .  (Lp  -  LL)  (T.L.)  =  99.11 

For  readings  29  3S  4,  5S  6  on  frame  one  -  above  procedure 

repeated. 

Thens       h±     =  99.11 

h,   =  99.17 

XS 

h.   =  99.13 

x3           K            =  99. 1^  =  h" 
h.   -  99.15      avg. 

x4 

h,   =  99.12 
5 

h.   =  99.13 
16 

Calculation  of  90^  confidence  interval  on  h~  of  frame  ones 

(h.   -  h")             (ht  -  h")2 
n                     n 

1                 -0.03                  0.0009 

2                 +0.03                  0.0009 

3                 -0.01                 0.0001 

4                +0.01                  0.0001 

5                -0.02                0.0004 

6                  +0.01                   0.0001 

2(b.   -  h")2=  0.0025 

' 
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Z(\     -   h")2 


S^  (h.) 


Z(\   -*") 


n  -  1 
2 


0.0025  =  0.0005 
5 


=  0.0005  =  0.000083 
6 


n(n-l) 

S    =  Ys2        =  Yo.  000083    =   0.0091 
0.05  =  Student's  t  for  959S  confidence  interval  =  2.571 
(t0.05)    (s)   =  2.571   (0.0091)   =  0.023^ 


).l4  +      0.02 


"0.05 


3.12  ^   E. 


99.16 


The  confidence  intervals  on  the  measurements  provided  a  means 
of  detecting  errors  that  occurred  at  any  time  during  the  processing 
of  the  data. 

Least  squares  line  (sample  regression  line): 


Observation  Number 


8 


•* 


99.14 


96 


03 


9\ 


192 


06 


2.\        2.9  2        Ze\ 


(Number  of  observations)  k  +  (Je)j  =£h. 

(2e)  k  +  (Ze2)  i    =£©h. 

Solve  simultaneous  equations  for  k  and  j. 
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Equation  of  least  squares  line: 
y  =  k  +  jx 

k  =  y  intercept  at  8  =  0 
j  =  slope  of  least  squares  line,  or 
sample  regression  coefficient 

Definition  of  Dynamic  Stability, Ifs     and Y f e* 

At  any  frame,  6,  the  absolute  deviation  from  regression  of  the 

height  of  a  disk  from  the  height  of  the  least  squares  line  through 

its  path  =  h.   -  (k  +  9j),  where 
19 

h   =  the  height  of  the  disk  at  6, 
*■• 
and  k  +  6j  ■  the  height  of  the  least  squares  line  at  8. 

Assumption:   Each  frame  over  the  descent  of  a  disk  was  equal 

to  the  same  increment  of  universal  time,  where 

the  conversion  factor  of  frames  to  seconds  was  F. 

For  any  randomly  selected  descent  of  a  disk  consisting  of  0 

frames,  the  summation  of  the  squares  of  absolute  deviation  from 

regression  = 

lh.   -  (k  +  j)    times  (number  of  frames  having  thia 

absolute  deviation,  which  was  always  one)  times 

(F)  +  |h±    -  (k  ♦  2j)|2   (1)  (F)  +  

....  +|hi2   -  (k+  »j)|2   (1)  (F). 
Units: 

|h.   -  (k  +  9j)|2  =  cm2. 

number  of  frames  having  above  absolute  deviation  =  frame. 

aeconds 
frames 
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2   . seconds 

cm.   **«•«.  t,^^ 

Therefore i 

Summation  of  absolute  deviations  from  regression  = 

\\    -  (k  ♦  ej)|2    (i)(f)  =  T&\&yJ. 

Finallyi 

fs2y<i.FZ|dyJ2/(e-2)(F) 

=  mean  square  deviation  from  regression,  or 

2 
dynamic  unstability  =   cm.   seconds 

frames   seconds 
frames 
•  cm. 

8  -  2  =  total  number  of  frames  minus  two,  since 
two  degrees  of  freedom  were  lost  in  the 
calculation. 

Dynamic  stability  =  U'  s     =        1 =       lr 

0r,  dynamic  stability  =  J  |s      =  ]W^\X       ■   -1 


Definition  of  Dynamic  Stability,  W 


AREA 


(See  PLATE  XIV,  Fig.  1,  where  the  shaded  area  on  this  plot  represents 

^f  Areae.) 
Trapezoidal  Rule: 
(See  PLATE  XIV,  Fig.  2  and  j) 

ab  =  |h±  -  (k  +  ©j)  |  On  Fig.  3« 

cd  -  I  h   -  fk  +  (9  +  l) jl  ab  +  cd  =  wy  -  iz 

le  +  'r  i  2 

ef  -  I  h.   -  [k  +  (8  +  2 )  jjl  Area  wyzx  -  area  abdc 

I  lax  a'1 


Le+J  ^.™..  seconds 


gh  =  |h   -  [k  +  (6  +  3)jJ|  yz  =  «i  =  (l)  F 

I   A*..*?  _      seconds 


seconds 


a 


Area  abdfhgeoa  =    Y  Area  =  ab  +  cd  (l)(F)  +  cd  +  ef  (l)(F) 

2  2 

♦  ef_±_gh  (1)(F)  =  (**  +  od  +  ef  +  §*)  *• 

Where  the  experimental  curve  crosses  the  least  squares  curve  be- 
tween consecutive  positions  of  the  disk  a  value  C  must  be  calculated. 
(See  PLATE  XIY,  Fig.  *f) 

{I- ,,-  c  .  ill'.  |S  J>.  <•  .  *%*} 

Area  a'b'c'd'a'   \ ^^  *  -  (< 

2{IVlk  +  9i,l   *  l1^.  7iLk  +  l9  +  1JdJ|3 

Situation  1.  (See   PLATE  XIII,    Fig.   2) 

Area  =  Area  al23it5678b 
fin.      -   (k  +  J)  I 
•\L1 -    +|h±     -  (k  +  2J)  |  +|h±     -  (k  +  3J)| 


♦  |h,    -(k  +  5j)|+  |hi6-(k  +  6i)|  +c6? 


5-  (k  +  7o)|         \K     -  2  (k  +  8j)| 


2 

Situation  2.  (See  PLATE  XIII,  Fig.  3) 

Area  =  Area  al23^5678b 


"I 


|\-(k  +  3)|     +  (s,  -  (k  ♦  u)j  +  |\-(k  +  33)l 

H-(k  +  ^)|2     |hi5-(k  +  ^)|   2 


2 


2 


h.   -  (k  +  73)       h   -  (k  +  8j) 

+  C,7  +  >_jj. 1  +   '_8 L 

•'  a  2 

Situation  3.   (See  PLATE  XIII,  Fig.  ^) 


+   c56  +  c6?  +  '_l 
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Area  =  Area  al23^56789b 


il 


\  '  ^  +  j)l      ♦   |hi2-(k+2,)|     +K ,-(k+3j) 


+     |h±     -  (k  +  4j)j    +|h.     -  (k  +  5d)|    +  \\     -  (k  +  6j)| 
Ih      -  (k  +  73)1    >.     -  (k  +  83)1        Ih      -  (k  ♦  93)U     N 

+      I    17  I     +l    18  '    +  J_^2 [(F), 


where   terms   'f  and  7   are   zero. 
In  summary: 

Area 


AREAe 


(cm.)    (fr^s)   If^nds 
(fra«e^) 
where  F  cancels   the  P  in(U,_E, 


(6-2)   P  (fra«e^)        seconds 

fraits» 


e 

l|'  AREA  =  Dynamic   stability  =     ),j^J}E 


Tp 


AREA 


Diagrams  of  Dynamic  Stability 
On  the  following  pages  diagrams  are  presented  illustrating  the 
methods  of  calculating  dynamic  stability. 


EXPLANATION  OF  PLATE  HII 

Definition  of  Dynamic  Instability. 
Fig.  1.  Physical  significance  of  Lg,  Lp,  and  Ll, 
on  single  frame  of  movie  of  disk  descent. 
Fig.  27| 


Three  possible  situations  in  calculation 
of  dynamic  ^instability,  J  AREA  . 


Fig.  A.J 
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PLiTE  XIII 
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Fig.  4 


EXPLANATION  OF  PLATE  XIV 

Definition  of  Dynamic  Instability,  IJJ  '  £&& 
Fig.  1.  Illustration  of^'iDEAg,  where  shaded  area  adjacent 

to  least  squares  line  represents Y'^gg^. 
Fig.  2.~)  Illustration  of  application  of  trapezoid  rule  to 
Fig.  3.J  calculating  area^'igg^. 
Fig.  U.     Illustration  of  calculation  of  area,  ^'jy^  >   ** 

special  case  where  trapezoid  rule  is  not  applicable. 
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PLATE  XIV 
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GENERAL  VISUAL  OBSERVATIONS  ON  THE  DESCENT 
OF  DISKS  THROUGH  TRIETHTLENE  GLYCOL  AND  WATER 


Discussion  of  Observations  and  Postulations 

Repeated  observations  were  made  of  the  general  behavior  of  the  disks 
as  each  descended  through  two  separate  fluid  media;  triethylene  glycol 
and  water.  For  water  seventeen  disks  were  investigated,  while  for 
triethylene  glycol  sixteen  were  investigated  due  to  the  buoyancy  of 
disk  one  in  the  glycol.  Preliminary  to  these  investigations  both  faces 
of  each  disk  were  marked  with  a  heavy  black  line  corresponding  to  the 
diameter  of  the  disk.  In  dropping  these  disks  through  either  triethylene 
glycol  or  water  the  lines  were  visible  to  an  observer  looking  down  through 
the  liquid  column  from  a  position  above  the  drop  tank.  From  this  position 
any  movement  of  the  diameter  of  the  disk  in  a  plane  parallel  to  the  sur- 
face of  the  earth  could  be  observed. 

The  disks  were  dropped  in  the  fluid  media  at  temperatures  of  80°  +  3°F 
for  water  and  90°  ±  3°F  for  triethylene  glycol.  At  80°F  the  density  of 
water  was  0.9963  and  at  90°F  that  of  one  hundred  per  cent  triethylene 
glycol  was  1.115.  However,  at  90°F  the  viscosity  of  triethylene  glycol 
was  22  centipoises  while  at  80°F  that  of  water  was  approximately  0.037 
times  this  or  0.82  centipoises. 

Observations  of  the  descents  of  each  disk  in  triethylene  glycol  showed 
that  the  marked  diameter  normally  did  not  rotate  throughout  the  descent. 
With  five  observations  of  each  disk  this  was  observed  in  96  per  cent  of 
all  the  observations.  Of  the  disks  in  triethylene  glycol  which  exhibited 
rotation  of  the  marked  diameter,  this  rotation  appeared  in  only  one  out  of 


48 


five  observations  made  on  each  of  these  disks.  In  all  cases  this 
rotation  was  observed  to  be  less  than  15°.  In  trie thy lene  glycol  the 
disk  occasionally  exhibited  rotational  motion  of  the  diametrical  axis  of 
maximum  slope,  but  even  in  the  majority  of  these  cases  the  marked  diameter 
remained  aligned  with  the  direction  which  it  assumed  when  released  (See 
Fig.  2,  PLATE  XVII).  These  observations  were  tabulated  in  Table  10, 
Appendix. 

In  water,  the  less  viscous  media,  the  marked  diameter  showed  de- 
cided deviations  from  the  direction  it  assumed  when  initially  released 
for  65  per  cent  of  the  one  hundred  observations.  In  water  at  80° ±  3°F 
two  of  the  heavier  disks,  30  and  32,  exhibited  no  diameter  rotation  over 
five  observations  of  their  descents.  The  remaining  fifteen  disks  ex- 
hibited rotation  of  the  marked  diameter  in  from  one  to  all  of  the 
observed  drops.  In  the  case  of  disks  19,  23,  4.0,   4-3,  and  45,  rotations 
first  in  one  direction  and  then  back  in  the  other  were  noted  in  part  of 
the  observations.  Also  it  was  observed,  by  the  utilization  of  a  second 
observer  watching  perpendicular  to  the  descent  of  the  disk,  that  many 
of  the  rotations  occurred  after  the  disk  had  descended  three-fourths  of 
the  depth  of  the  column.  These  observations  were  tabulated  in  Table  11, 
Appendix. 

Other  observations  were  made  on  the  path  of  motion  which  the  center 
of  the  disk  followed  as  it  descended.  When  viewed  from  above  the  liquid 
column  these  obsprvations  showed  that  in  ninety-seven  per  cent  of  five 
descents  each  of  sixteen  different  disks  the  path  of  a  given  disk  re- 
mained within  one  vertical  plane.  A  sighting  mechanism,  built  of  welding 
bars,  was  constructed  in  the  shape  of  a  circle,  the  same  size  as  the  top 
of  the  drop  column,  and  with  a  single  bar  as  a  diameter  of  the  circle. 
The  apparatus  was  fitted  on  the  top  of  the  drop  column  such  that  as  a 
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disk  descended  it  was  possible  to  rotate  the  mechanism  so  that  the  dia- 
metrical bar  was  contained  in  the  vertical  plane  through  the  path  of 
motion  of  the  descending  disk.  A  pointer  was  added  to  the  mechanism 
such  that  it  pointed  to  values  of  angle  inscribed  along  the  rim  of  the 
drop  column.  The  angle  or  angles  between  the  vertical  plane  containing 
the  disk  path  of  motion  and  a  vertical  reference  plane  passing  through 
0°  on  the  rim  of  the  column  were  determined  for  four  separate  drops  in 
triethylene  glycol.  These  observations  in  Table  12,  Appendix,  showed 
that  only  2  descents  out  of  64  exhibited  a  shift  in  the  vertical  plane 
passing  through  the  path  of  motion  of  the  disk  through  its  100  cm. 
descent. 

No  such  conclusive  observations  were  made  for  the  seventeen  disks 
descending  in  water.  In  water  57  per  cent  of  the  descents  showed  no 
single  planar  restriction  of  the  path  of  motion  as  observed  from  the  top 
of  the  column.  Of  the  43  per  cent  of  the  descents  showing  single  planar 
motion  73  per  cent  of  these  observations  were  made  on  disks  20,  23,  27, 
30,  32,  35,  37,  and  40.  In  32  per  cent  of  the  single  planar  motions  the 
disk  hit  the  side  of  the  tank  before  the  full  100  cm.  between  reference 
lines  had  been  traveled.  In  no  case  did  any  disk  show  single  planar 
wave,  flutter,  or  oscillatory  motion  for  all  five  observations  in  water 
at  80° ±  3°F.  These  observations  were  tabulated  in  Table  14,  Appendix. 

The  observations  in  this  study  were  restricted  to  disks  descending 
in  triethylene  glycol  and  in  water  at  Reynolds  numbers  above  78.  After 
a  multitude  of  observations,  including  those  presented  in  Tables  13  and 
14  of  the  Appendix,  it  was  postulated  that  the  disks  descended  in  one  or 
more  of  four  main  types  of  motion.  These  motions  were  named  stable, 
flutter,  wave,  and  spiral. 
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In  stable  motion  the  center  of  the  disk  fell  along  a  vertical  line. 
The  vertical  line  was  always  perpendicular  to  the  face  of  the  disk.  In 
triethylene  glycol  particles  6  and  46  descended  in  a  stable  manner  in 
five  out  of  five  descents  of  each  disk.  In  water  no  disk  exhibited 
stable  descent  in  any  of  the  five  drops  tabulated  in  Table  14,  Appendix. 
At  least  ten  other  drops  not  tabulated  substantiated  these  observations 
over  the  same  range  of  disks. 

In  flutter  the  disk  exhibited  what  appeared  to  be  a  wobble  or 
"teeter-totter"  motion  about  a  given  fulcrum.  This  fulcrum  consisted 
of  one  particular  diameter  of  the  disk;  this  diameter  did  not  change  on 
the  disk  throughout  the  descent.  Flutter  occurred  while  the  center  of 
the  disk  was  descending  through  either  a  vertical  line  or  a  three 
dimensional  curvelinear  path.  In  triethylene  glycol  this  diametrical 
fulcrum  remained  constant  in  direction  throughout  the  descent. 

In  triethylene  glycol  disks  9,  13,  20,  and  45  exhibited  what  was 
observed  as  flutter  in  one,  two,  or  three  descents  out  of  five.  Disks 
19,  20,  and  30  exhibited  flutter  as  a  component  of  wave  motion.  In 
water  disks  13,  16,  19  and  20  appeared  to  descend  with  a  flutter  in  at 
least  one  observation  out  of  seven.  In  Table  14,  Appendix,  only  five  of 
these  observations  were  tabulated.  It  was  found  in  the  case  of  flutter 
that  a  second  or  third  observer  sometimes  viewed  the  motion  as  a  small 
wave,  instead  of  flutter.  This  led  to  later  decisions  concerning  flutter 
as  a  postulated  type  of  motion. 

Wave  motion  was  described  as  an  oscillation  of  the  center  of  the 
disk  with  the  face  of  the  disk  always  remaining  perpendicular  to  the 
direction  of  motion  (See  Fig.  3,  PLATE  XV).  In  oscillation  the  path 
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of  motion  of  the  center  of  the  particle,  when  projected  onto  a  vertical 
plane,  appeared  similar  to  a  sine  or  cosine  curve.  As  already  described 
for  descents  in  triethylene  glycol  a  single  vertical  plane  contained 
the  entire  path  of  oscillatory  motion.  For  water  this  restriction  to  a 
single  vertical  plane  did  not  apply,  as  oscillation  and  other  motions  in 
water  were  unmistakably  three  dimensional,,  In  triethylene  glycol  75  per 
cent  of  the  80  observations  in  Table  13,  Appendix,  revealed  wave  motion 
alone  or  in  combination  with  flutter  or  disk  rotation  without  diameter 
rotation.  In  water  wave  motion  alone,  or  in  combination  with  spiral  or 
disk  rotation  with  diameter  rotation,  appeared  in  67  per  cent  of  the  85 
observations  of  3  drops  in  Table  14,  Appendix.  Spiral  was  defined  as  a 
"cork  screw"  path  of  motion  of  the  center  of  the  disk  (See  Fig.  1, 
PLATE  XVII).  In  Table  13,  Appendix,  spiral  motion  did  not  appear  in  5 
drops  through  triethylene  glycol.  Other  observations  made  in  this  study, 
but  not  tabulated,  substantiated  this.  In  water  disk  one  exhibited  ex- 
clusively spiral  motion  over  three  drops  as  shown  in  Table  14,  Appendix. 
In  29  per  cent  of  the  total  observations  of  seventeen  disks  in  Table  14, 
Appendix,  spiral  motion  appeared  either  alone  or  as  an  extension  of  wave 
motion  and  flutter.  The  two  dimensional  appearances  of  these  motions 
were  illustrated  in  PLATES  XV  and  XVII . 

Other  postulated  types  of  motion  were  also  recognized  but  these 
motions  only  seldom  or  never  appeared  in  triethylene  glycol.  The  nomen- 
clature given  two  of  these  motions  were  disk  rotation  without  diameter 
rotation  and  disk  rotation  with  diameter  rotation.  In  both  cases  disk 
rotation  was  basically  defined  as  rotation  of  the  diametrical  axis  of 
maximum  slope  of  a  disk.  In  the  former  case  the  disk  did  not  actually 
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rotate,  but,  without  a  diameter  reference  mark,  appeared  to  do  so  (See 
Fig.  2,  PLATE  XVII).  In  the  latter  case  the  disk  rotated  as  apparent  from 
the  rotation  of  a  marked  diameter  on  the  face  of  the  disk  (See  Fig.  3, 
PLATE  XVII),  It  was  found  that  some  discernment  was  necessary  in  dis- 
tinguishing between  borderline  cases  of  spiral,  disk  rotation,  and  wave. 

If  the  motion,  appearing  as  a  wave  over  a  small  increment  of  the 
total  descent,  exhibited  at  the  same  time  one  complete  loop  of  a  spiral 
over  the  total  descent  it  was  then  classified  as  spiral  instead  of  wave. 
If  the  motion  appearing  as  a  disk  rotation  exhibited  any  recognizable 
segment  of  a  spiral  loop  over  the  total  descent,  it  was  classified  as 
a  spiral.  In  the  ease  of  disk  rotation  and  wave  motion  together,  the 
combination  was  treated  as  wave  motion  with  rotation.  Infrequently,  in 
only  5  per  cent  of  the  observations  in  Table  13,  Appendix,  disk  rotation 
without  diameter  rotation  appeared  in  triethylene  glycol,  while  disk 
rotation  with  diameter  rotation  did  not  appear  in  any  of  the  80  ob- 
servations. In  water,  disk  rotation  with  diameter  rotation  appeared 
in  less  than  5  per  cent  of  the  observations  in  Table  14,  Appendix,  while 
disk  rotation  without  diameter  rotation  was  not  apparent  in  any  of  the 
85  observations.  Diameter  rotation  in  water  was  most  closely  associated 
with  spiral  motion  as  was  seen  from  a  comparison  of  Tables  11  and  H, 
Appendix.  These  were  illustrated  in  PLATE  XVII. 

At  this  point  it  was  discovered  that  fairly  accurate  reproductions 
of  the  disk  path  of  descent  could  be  obtained.  This  was  accomplished 
by  tracing  with  a  pencil  the  projected  image  of  the  disk  from  each 
movie  frame,  as  each  frame  in  turn,  was  projected  onto  a  screen  of 
white  paper.  Prior  to  tracing,  each  frame  was  adjusted  to  the  position 
on  the  screen  which  corresponded  to  that  previously  occupied  by  the 
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preceding  frame  traced.  This  positioning  was  accomplished  by  adjusting 
the  upper  reference  line  pictured  on  each  frame  of  the  movie  to  a  con- 
stant reference  point  on  the  white  paper  background;  by  this  technique 
a  fairly  accurate  reproduction  of  the  overall  path  of  descent  was  ob- 
tained. Only  parts  of  the  six  descents  filmed  were  analyzed  in  this 
manner  due  to  the  excess  time  and  personnel  requirements. 

From  these  sketches  more  exacting  pictures  of  the  postulated  types 
of  motions  of  disks  in  triethylene  glycol  were  established.  The  results 
of  analysis  of  these  sketches  were  tabulated  in  Table  15,  Appendix.  It 
was  found  that  flutter  was  merely  wave  motion  with  barely  perceptible 
horizontal  displacement  when  viewed  at  the  true  velocity  of  descent. 
This  observation  was  further  substantiated  by  analyzing  the  remainder 
of  the  movies.  For  this  analysis  a  vertical  reference  line  was  estab- 
lished by  holding  a  white  string  vertically  downward  from  the  initial 
projected  position  of  a  disk  as  viewed  on  the  movie  screen.  The  pro- 
jection of  the  disk  was  carefully  observed  as  it  descended  along  the 
string  on  the  screen.  Horizontal  displacements  to  the  right  and  left 
of  this  line  Immediately  showed  up  on  all  the  disks  which  had  been 
formerly  depicted  as  descending  with  a  fluttering  motion.  These  same 
operations,  also,  showed  that  the  disks  depicted  as  settling  in  stable 
motion  did  indeed  settle  in  that  manner.  Typical  results  of  these 
latter  analyses  were  also  tabulated  in  Table  15,  Appendix. 

From  these  newer  approaches  to  observing  motion  it  was  postulated 
that  wave  motion  actually  consisted  of  three  types.  That  previously 
classified  as  wave  motion  was  reclassified  as  two  types  of  motion  - 
wave  I  and  wave  III.  Wave  I  was  subtyped  as  wave  I—A  and  wave  I-B. 
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Wave,  as  previously  defined,  became  wave  I-A.  It  was  noted  that  wave 
I-A  consisted  of  planar  oscillatory  motion  and  it  was  postulated  that 
the  disk  always  maintained  orientation  such  as  to  offer  the  face  of 
greatest  resistance  perpendicular  to  the  path  of  motion  of  the  center 
of  the  disk.  Wave  I-B  motion  was  observed  as  similar  in  definition  to 
wave  I-A  motion,  but  without  the  restriction  that  the  face  of  greatest 
resistance  always  be  perpendicular  to  the  path  of  the  center  of  the 
disk.  Analysis  of  sketches  of  disk  descents  showed  that  as  a  disk 
became  heavier  its  increased  moment  of  inertia  prevented  it  from  be- 
having as  a  wave  I-A  disk.  It  was  impossible  to  discern  between  these 
two  types  of  motion  on  the  sketches  of  the  disk  descent.  However,  as  the 
moment  of  inerita  of  the  disks  increased,  wave  I-A  and  wave  I-B  motions 
formed  a  logical  sequence  between  stable  and  wave  III  motions.  Wave  I-B 
motion  wa3  definitely  observed  in  disks  19,  20,  23,  and  27.  However, 
for  disks  9,  13,  and  16  there  was  little  to  discern  between  the  types  of 
motion,  since  the  sketches  were  not  refined  enough  to  allow  precision 
study.  The  heavier  a  disk  became  the  larger  the  absolute  difference  of 
9  -  90°  became.  8  represented,  at  any  instant,  the  smallest  angle  between 
the  plane  of  the  leading  face  of  the  disk  and  the  directional  vector  of 
the  center  of  the  disk.  Disks  exhibiting  wave  I-A  and  I-B  motion,  at 
least  once,  were  16,  19,  20,  23,  27,  and  30.  These  motions  are  illustrated 
in  PLATES  XV,  XVI,  and  XVII. 

When  the  plane  containing  the  face  of  a  disk  rotated  away  from  the 
horizontal  through  the  vertical  and  back  to  the  horizontal  wave  III  motion 
was  observed.  It  was  postulated  that  the  angle  between  the  plane  con- 
taining the  face  of  the  disk  and  a  horizontal  plane  became  increasingly 
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larger  as  the  disk  descended  until  the  disk  turned  over  or  flipped.  Flip 
was  defined  as  the  motion  of  a  disk  where  the  lower  face  of  the  disk  inter- 
changed with  the  upper  face  through  a  180°  rotation.  Disks  exhibiting 
wave  III  motion  were  disks  with  large  moments  of  inertia  with  respect  to 
their  weights.  In  all  cases  they  were  of  steel  rim  and  plastic  core 
construction.  Therefore  this  flipping  action  was  a  result  of  the  large 
moments  of  inertia  of  these  disks.  This  was  a  reasonable  conclusion  when 
it  was  considered  that  a  baton  is  constructed  of  a  shaft  with  heavy  weights 
at  each  end  such  that  it  has  a  large  moment  of  inertia. 

The  motions  described  here  and  illustrated  in  PLATE  XVI  can  be  ob- 
served by  the  use  of  a  sheet  of  typing  paper.  A  rectangle  U  by  l£  inches 
cut  from  the  sheet  and  dropped  face  down  from  shoulder  high  exhibits  wave 
III  motion.  A  square  consisting  of  three  2£  by  2J-  pieces  glued  together 
exhibits  either  wave  I-A  or  wave  I-B  motions. 

Further  observations  concerned  the  frequencies  of  planar  oscillation 
of  disks  descending  in  wave  patterns.  From  the  sketches  prepared  for 
Table  15,  Appendix,  the  following  frequencies  in  cycles  per  100  cm. 
descent  were  observed  for  parts  of  two  different  drops. 

Table  Is  Frequencies  in  Cycles  per 
100  cm.  Descents  of  Disks. 


Disk  Number 

Drop  One 

Droo  Two 

9 



6 

13 



6 

16 

6 

6 

19 

6 

5* 

20 



5  3/5 

23 

5 



27 

4  3/4 

. 

30 

U 



32 

U 



35 

_ _ 

__- 

37 

3  2/3 

■ — ■ 

40 
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Diagrams  of  Postulated  Types  of  Motion 


On  the  following  pages  diagrams  are  presented  illustrating  the 
postulated  types  of  motion  of  disks  descending  in  triethylene  glycol 
and  water. 


EXPLANATION  OF  PLATE  XVI 

Postulated  types  of  motion  of  disks  descending  in  a  fluid  media. 
Fig.  1.  Wave  I-E,  the  other  subtypes  of  Wave  I  motion. 
Fig.  2.  Wave  III. 
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PLATE  XVI 


WAVE      l-B 


WAVE 
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EVALUATION  OF  DATA  AND  PRESENTATION  OF  RESULTS 

Dynamic  stabilities  were  calculated  by  three  methods;  each  method 
was  previously  described.  The  dynamic  stabilities  evaluated  were 
^s2_  x,  reciprocal  of  mean  square  deviation  from  regression,  t1s  y„x> 
reciprocal  of  sample  standard  deviation  from  regression,  and  IjJ  AREA? 
reciprocal  of  mean  area  between  least  squares  line  and  path  of  disk 
descent.  These  values  for  three  drops  were  tabulated  in  Tables  7,  8, 
and  9  of  the  Appendix.  Obvious  differences  existed  between  the  dynamic 
stabilities  calculated  by  any  one  method  for  each  disk  over  its  three 
descents.  The  more  outstanding  reasons  for  these  differences  were  the 
randomness  of  free-fall  descents  of  disks  in  fluids,  the  paralactic 
errors,  the  measuring  errors,  and  the  extraneous  disturbances,  as 
described  in  the  procedure,  in  the  fluid  media  during  descents  of  the 
disks. 

It  was  impossible  to  distinguish  between  true  randomness  of  motion 
and  erraticness  of  motion  which  was  induced  by  extraneous  disturbances. 
Extraneously  induced  disturbances  in  the  fluid  were  assumed  to  be 
sufficiently  dampened  by  the  precautions  earlier  described  in  the  pro- 
cedure, although  these  steps  certainly  guaranteed  no  possibility  for 
complete  removal  of  this  source  of  error.  Measuring  errors,  a  result 
of  obtaining  values  of  disk  height  from  motion  pictures  of  disk  descents, 
provided  measureable  but  unpredictable  sources  of  disagreement,  dictated 
by  equipment,  personal,  and  environmental  deficiencies  at  any  time. 
Parallax  -  as  treated  by  Himes  (3)  -  also  afforded  measureable  sources  of 
error,  where  the  magnitude  of  paralactic  error  was  determined,  for  the 


EXPLANATION  OF  PLATE  XIX 

Illustration  of  90?  confidence  intervals 

on  *f  I  s  '   ,  dynamic  stability,  for  range 

of  disks  investigated. 

""A.   upper  and  lower 
v  limits  of  90? 
1  confidence  Intervals. 

X   Drop  1. 

O   Drop  2. 

D   Drop  3. 
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PLATE  XIX 
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EXPLANATION  OF  PLATE  XX 

Plots  of  it3  y  x»  dynamic  stability,  for  three 
drops  over  entire  range  of  disks  showing  inten- 
sities of  y^s   x  in  relation  to  postulated 
types  of  motion. 
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EXPLANATION  OF  PLATE  XXI 

Plots  of  (|/  ajjeas  dynamic  stability,  for  three 
drops  over  entire  range  of  disks. 
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EXPLANATION  OF  PLATE  XHI 

Plots  of  Ys2y.x?  dynamic  stabilities  for 
three  drops  over  entire  range  of  disks. 
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EXPLANATION  OF  PLATE  XHII 

Bar  graph  illustrating  relationship  of  9056 
confidence  intervals  of  ffs^   ,  dynamic 
stabilities,  of  disks  classified  according 
to  postulated  types  of  motion. 
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major  part,  by  the  position  of  the  particle  relative  to  the  location  of 
the  reference  wires  and  to  the  location  of  the  camera.  By  moving  the 
camera  to  a  position  further  from  the  drop  column  the  parallax  was 
lessened  by  a  calculable  amount.  However,  the  allowable  distance  be- 
tween the  camera  and  the  tank  was  limited  by  the  characteristics  of  the 
available  telescopic  lenses. 

In  summary  extraneously  induced  disturbances  of  the  descent  of  a 
disk  and  randomness  of  behavior  in  the  descent  of  a  disk  provided 
incalculable  and  uncontrollable  sources  of  disagreement  between 
measured  paths  of  descents.  Measuring  error  was  calculable  but  un- 
controllable, while  parallax  was  both  calculable  and  controllable  within 
equipment  and  camera  limitations.  Therefore,  it  remained  to  minimize 
paralactic  error  to  the  degree  where  it  was  of  the  same  or  of  a  lesser 
magnitude  than  measuring  error.  Without  altering  the  original  locations 
of  the  reference  wires  this  was  accomplished  by  increasing  the  distance 
between  the  camera  and  the  drop  column  to  the  maximum  allowed  by  the 
available  telescopic  lenses.  Confidence  intervals  on  maximum  measuring 
error  when  compared  with  maximum  possible  paralactic  error  were  found 
to  be  of  the  same  range  of  magnitude.  This  comparison  was  presented  in 
PARALLAX  (3)  in  conjunction  with  a  rigorous  investigation  of  parallax. 

Other  sources  of  error  were  due  to  the  non-uniformity  of  release 
of  the  disks  and  size  limitations  of  the  drop  column.  The  former  was 
minimized  by  disregarding  an  initial  length  of  each  descent  while  the 
latter  was  unavoidable  for  monetary  reasons.  Another  source  of  error 
resulted  from  variations  in  the  rate  at  which  the  camera  filmed  the 
descent.  The  speed  of  the  camera  was  found  to  vary  from  16.5  to  19.6 
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frames  per  second.  This  meant  that  the  relative  magnitude  of  each 
frame  varied  within  a  given  descent.  The  assumption  that  all  frames 
were  equal  to  the  same  increment  of  time  was  not  true.  This  variation 
was  a  characteristic  of  the  spring  and  governor  regulating  the  camera 
speed.  It,  also,  depended  upon  how  often  and  to  what  degree  the 
photographer  wound  the  camera.  For  this  reason  it  was  impossible  to 
calculate  the  magnitude  of  this  source  of  error.  The  solution  to  this 
problem  lay  in  the  statement,  "An  ounce  of  prevention  is  worth  a  pound 
of  cure".  However,  "the  ounce  of  prevention,"  consisted  of  a  synchro- 
nizing motor  attachment  for  the  camera.  Monetary  limitations  prevented 
this.  However,  after  some  consideration  it  was  felt  that  this  lack  of 
uniformity  in  the  camera  speed  was  not  nearly  as  critical  as  measuring 
and  paralactic  errors,  since  the  variation  of  camera  speed  within  any 
one  descent  was  extremely  small.  This  wa3  borne  out  by  the  excellent 
agreement  between  the  calculated  dynamic  stabilities  of  various  descents 
of  most  disks. 

The  values  of  dynamic  stability  by  each  method  -  Y*s  y.x> 
f  AREA'  y  s2J.x>  ~   were  Plotted  in  PLATES  XX,  XXI,  and  XXII,  re- 
spectively. The  method  based  upon  Tf  s  y.x  was  randomly  chosen  from 
these  three  for  further  consideration.  In  PLATE  XIX,  ^ferlTx  was 
plotted  as  the  upper  and  lower  limits  of  90  per  cent  confidence  intervals 
oil;  the  dynamic  stabilities  of  each  disk.  This  graph  illustrated  the 
statement  that  the  true  value  of  dynamic  stability  for  each  disk  had  a 
90  per  cent  chance  of  lying  between  the  individual  limits  described. 
This  did  not  state  that  the  true  curve  of  dynamic  stability  had  a  90 
per  cent  chance  of  lying  entirely  between  these  limits  formed  by 
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connecting  the  individual  limits.  However,  the  limits  so  acquired  by 
connecting  these  individual  limits  lead  to  the  statement  that  each 
point  of  the  true  curve  of  dynamic  stabilities  had  a  90  per  cent  chance 
of  lying  within  these  limits. 

During  the  period  of  calculation  of  dynamic  stabilities  a  completely 
isolated  study  was  being  conducted  on  the  types  of  motion  apparent  in  the 
disk  descents.  As  was  described  earlier  in  this  thesis  under  General  Visual 
Observations  on  the  Descent  of  Disks  through  Triethylene  Glycol  and 
Water  each  disk  was  assigned  to  one  of  three  main  types  of  motion.  In 
PLATE  XX  the  ranges  of  assigned  motions  were  compared  with  the  magnitudes 
of  dynamic  stabilities.  In  PLATE  XXIV,  a  similar  comparison  was  made  be- 
tween physical  characteristics  of  the  disks  and  assigned  motions.  In 
PLATE  XXV,  a  like  comparison  was  made  between  Reynolds  numbers  and 
assigned  motions. 

DISCUSSION  OF  RESULTS 

With  reference  to  Tables  7,  8,  and  9  of  the  Appendix,  it  is  seen 
that  decided  agreements  between  dynamic  stabilities  of  the  three  descents 
of  each  disk  were  clearly  observed  for  most  of  the  disks.  Disks  not 
exhibiting  this  agreement  were  numbers  9,  30,  32,  and  45.  In  these  cases 
the  disagreements,  relative  to  the  magnitudes  of  the  dynamic  stabilities, 
were  large.  Erratic  behavior  in  the  source  of  a  measurable  quantity  was 
associated  by  this  investigator  with  the  inability  to  reproduce  measure- 
ments of  that  quantity.  Therefore,  because  of  the  visually  observed 
erraticalness  of  disks  30,  32,  35,  37,  and  40  during  their  descents  dis- 
agreement was  expected,  although  only  disks  30  and  32  actually  exhibited 
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any  marked  diversion.  This  indicated  some  degree  of  uniformity  of  be- 
havior within  the  erraticalness.  Erraticalness  in  the  descent  of  a  disk 
was  also  associated  with  low  dynamic  stabilities.  This  was  true;  disks 
32,  35,  27,  and  40  had  the  lowest  values  of  dynamic  stability. 

Dynamic  stabilities,  ^s2   ,  ^  AREA,  and  *f  s2y#x,  were  plotted  on 
PLATES  XX,  XXI,  and  XXII,  respectively.  From  the  plots  of  all  three 
methods  four  levels  of  magnitudes  of  dynamic  stability  were  apparent. 
Selecting  f^s^    for  further  study,  it  was  found  that  disk  6,  the  most 
stable,  formed  the  highest  level.  Disks  13,  16,  and  43  formed  the  second 
level.  Disks  19,  20,  23,  and  27  formed  the  third.  Disks  32,  25,  27,  and 
40  formed  the  lowest  and  final  level.  Disks  9,  30,  and  45  had  values  of 
dynamic  stability  in  more  than  one  level.  With  only  two  values  of  dynamic 
stability  it  was  improper  to  classify  disk  46  in  any  of  the  levels. 

The  four  levels  were  then  compared  with  the  postulated  types  of 
motion  exhibited  by  the  descending  disks,  Table  15,  Appendix.  This  study, 
dealing  with  the  types  of  motion,  was  conducted  prior  to  and  separate 
from  the  evaluation  of  dynamic  stabilities.  Therefore,  the  excellent 
agreement  of  the  two  was  surprising.  These  postulated  types  of  motion 
were  compared  with  dynamic  stabilities  on  the  plot  of  WVs*        against 
disk  number,  PLATE  XX.  The  discrepancies  caused  by  disks  9,  30,  and 
45  were  resolved  by  considering  these  disks  as  transitional  between 
two  types  of  motion.  Table  15,  Appendix,  substantiated  this  for  disk 
30,  since  this  disk  was  found  to  exhibit  more  than  one  type  of  motion. 
Disk  43  was  observed  to  exhibit  wave  I.  Its  values  of  dynamic  stability 
did  not,  however,  indicate  to  which  subtype  it  belonged.  Since  one  value 
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of  dynamic  stability  on  disk  43  fell  into  the  wave  I-A  level  and  the  other 
two  into  the  wave  I-B  level,  disk  43  was  finally  classified  as  transitional 
between  wave  I-A  and  wave  I-B.  Disk  9  presented  the  only  outstanding  dis- 
crepancy in  that  the  observed  motion,  postulated  as  wave  I,  did  not  agree 
with  the  levels  in  which  its  dynamic  stabilities  fell  -  wave  I-B  and  wave 
III.  Therefore,  disk  9  was  classified  as  being  in  an  undeterminable  trans- 
ition region. 

Ninety  per  cent  confidence  intervals  were  calculated  on  the  three 
values  of  dynamic  stability  for  each  disk,  as  presented  in  Table  3, 
Appendix.  The  close  agreement  between  the  confidence  intervals  on 
dynamic  stability  among  various  disks  within  a  given  type  of  fall  in- 
dicated that  the  dynamic  stabilities  for  these  disks  came  from  the  same 
populations.  The  disks  having  comparable  confidence  intervals  were 
found  to  form  four  separate  groups.  These  four  groups  corresponded 
exactly  to  the  four  levels  of  dynamic  stability  earlier  found  to 
correspond  to  the  postulated  types  of  motion.  Therefore,  the  values  of 
dynamic  stability  for  each  disk,  with  the  exceptions  of  9,  30,  and  45> 
fell  into  one  of  the  postulated  types  of  motion.  Ninety  per  cent  con- 
fidence intervals  were  then  calculated  on  each  group  as  presented  in  Table 
4,  Appendix.  The  limits  on  these  groups,  or  levels,  were  then  presented 
as  a  bar  graph  in  PLATE  XXIII.  This  illustrated  the  statement  that  the 
true  value  of  dynamic  stability  for  each  of  these  groups  had  a  90  per 
cent  chance  of  lying  within  the  limits,  indicating  that  each  such  level 
was  distinct  and  separate  from  the  others. 

The  next  step  was  to  make  an  attempt  to  correlate  the  properties  of 
the  disks  and  fluid  to  the  assigned  types  of  motion.  The  properties  of 
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the  disks  were  tabulated  in  Table  5,  Appendix,  and  plotted  in  conjuction 
with  the  types  of  motion  assigned  each  disk  in  PLATE  XXIV.  From  Tables 
7,  8,  and  9  of  the  Appendix  Re's  were  plotted  in  the  same  manner  in 
PLATE  XXV.  From  PLATE  XXIV  the  only  conclusive  statement  possible  was 
that  disks  exhibited  greater  instability  as  their  moments  of  inertia 
became  large  relative  to  their  weights.  Stable  motion  occurred  at  low 
weights,  low  moments  of  inertia,  and  low  densities  relative  to  the  com- 
parable properties  for  disks  which  exhibited  other  postulated  types  of 
motion.  Disk  6  exhibited  stable  motion.  Finally,  from  PLATE  XXV,  the 
Reynolds  numbers,  based  on  Ds,  for  the  three  descents  of  disk  6  were  less 
than  100.  Disks  exhibiting  waves  I-A,  I-B,  and  III  motion  had  Re's  of 
275  to  325,  225  to  650,  and  300  to  625,  respectively. 

CONCLUSIONS 

The  most  significant  conclusion,  based  upon  more  than  41,000  measure- 
ments, was  that  dynamic  stability  of  disks  is  measureable  and  therefore 
exists. 

It  was  concluded  that  a  disk,  in  a  free-fall  descent,  exhibited  one 
of  three  postulated  types  of  motion;  i.  e.  stable,  wave  I,  or  wave  III. 
It  was  further  postulated  that  wave  I  motion  consisted  of  subtypes;  wave 
I-A  and  wave  I-B.  The  distinction  between  these  subtypes  appeared  in 
the  measured  dynamic  stabilities;  corresponding  to  the  exact  boundary 
predicted  by  the  observations  (See  Table  15,  Appendix).  In  water  disks 
exhibited  three  main  types  of  motion;  i.  e.  wave,  spiral,  or  disk  rotation 
with  diameter  rotation,  or  combinations  of  these.  For  triethylene  glycol 
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a  definite  relationship  was  found  between  these  motions  and  dynamic 
stability,  lffs2y#x. 

Table  2:  Relationship  Between  Types  of  Motion  and  ^l3  y.x« 


Postulated 

:   Range  of  dynamic 

Disk  number     s 

type  of  motion 

:   stabilities,  l/cm. 

6 

Stable 

11.49  to  14.92 

Tl, 

16,  and  43 

Wave  I-A 

6.02  to  8.55 

1q, 

20,  23,  and  27 

Wave  I-B 

4.15  to  5.46 

32, 

35,  37  and  40 

Wave  III 

0.44  to  2.01 

Each  range  of  dynamic  stabilities  corresponding  to  a  particular  type  or 
subtype  of  motion  was  distinct  and  separate  from  any  of  the  others.  By 
applying  confidence  intervals  to  these  ranges  it  was  shown  that  the  true 
value  of  dynamic  stability  for  each  type  of  motion  had  a  90  per  cent 
chance  of  lying  within  limits  which  did  not  overlap  the  limits  of  other 
types  of  motion. 

Actually,  three  methods,  IWE2   ,p7I,  and  ty  AREA  were  used  in 
calculating  dynamic  stability.  The  method  based  upon  (Ms 2    was  randomly 
chosen  from  these  three  for  further  consideration. 

From  a  comparison  of  physical  properties  of  the  disks  against  the 
types  of  motion  into  which  each  disk  fell  it  was  found  that  a  low  dynamic 
stability  was  associated  with  a  high  moment  of  inertia  relative  to  the 
weight  of  the  disk.  This  is  summarized  in  PLATE  XXIV. 

From  a  comparison  of  Reynolds  numbers  for  each  disk  with  the  types 
of  motion  of  each  disk  it  was  concluded  that  a  Re,  based  on  Dg,   of  one 
hundred  or  less  was  associated  with  stable  motion.  Disk  6,  the  only 
disk  which  clearly  exhibited  stable  motion,  was  the  lightest  one  and  it 
also  had  the  least  density  and  moment  of  inertia  of  all  those  investigated. 
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These  conclusions  were  all  restricted  to  disks  which  descended  in 
triethylene  glycol  at  88  to  90°F  at  Re,  based  on  Ds,  ranging  from  78 
to  640. 

RECOMMENDATIONS 

The  following  are  recommended  for  immediate  investigations 

1.  Paralactic  error  may  be  noticeably  reduced  by  moving  the  refer- 
ence wires  to  a  position  in  contact  with  the  drop  column.  This  suggests 
an  investigation  of  the  positioning  of  the  top  reference  wire  immediately 
above  the  release  mechanism  and  the  bottom  reference  wire  inside  the  column 
near  the  bottom. 

2.  Parallax  may  be  further  reduced  by  removing  the  camera  further 
from  the  drop  column,  thus  suggesting  that  an  investigation  be  made  of 
telescopic  lenses  allowing  this  increased  distance. 

3.  Recent  films  of  disk  descents  revealed  that  a  suspension  of  rust 
pigments  in  the  triethylene  glycol  produced  more  distinct  pictures  of  the 
disks.  Other  pigments  may  prove  even  more  satisfactory. 

U.     Error,  although  not  significant  when  compared  to  paralactic  and 
measuring  errors,  was  introduced  into  the  final  values  of  dynamic 
stability  due  to  the  non-uniformity  of  the  camera  speed.  It  is  possible 
to  rectify  this  by  the  use  of  a  synchronizing  motor.  As  paralactic  and 
measuring  errors  are  decreased  this  consideration  will  become  important. 

5.  Measuring  errors  and  measuring  time  could  be  noticeably  de- 
creased by  moving  the  dark  room  and  projector  to  a  location  supported 
by  concrete  flooring. 


6.  The  size  of  the  drop  column  should  be  appreciably  increased.  An 
increased  diameter  would  lessen  wall  effects,  while  both  an  increased 
diameter  and  length  would  allow  a  more  thorough  study  of  disks  exhibiting 
wave  III  motion. 

The  following  are  recommended  for  long-range  investigation' 

1.  An  analysis  should  be  made  of  the  possibility  of  using  a  radio- 
active source  to  expose  film  strips  attached  to  the  drop  column  to  obtain 
a  continuous  picture  of  the  disk  descent. 

2.  A  survey  should  be  made  of  the  possibility  of  employing  either 
electric  eyes  or  magnetic  coils  along  the  length  of  the  column,  such  that 
with  the  addition  of  an  appropriate  recorder,  readings  of  particle  height 
may  be  directly  obtained. 

3.  A  study  should  be  made  of  drop  tanks  other  than  cylindrical  in 
shape. 

Statements  could  be  made  concerning  what  next  should  be  investigated, 
but  dynamic  stability  has  revealed  such  a  fertile  field  as  to  make  any 
recommendations  for  future  work  little  more  than  twigs  in  an  ocean, 
therefore,  such  concern  will  be  left  to  the  interests  of  future  investi- 
gators . 
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Triethylene  Glycol  as  the  Fluid  Media 

Triethylene  Glycol  was  a  relatively  new  compound  commercially  at  the 
initiation  of  this  research.  Its  high  boiling  point,  low  volatility,  low 
toxicity,  low  corrosion  rate,  and  high  stability  made  it  ideal  as  a  high 
boiling  solvent.  These  same  properties  plus  its  high  viscosity  at  room 
temperatures  and  its  transparency  made  it  similarly  ideal  for  the  study 
of  particle  motions.  These  properties  are  presented  in  tables  and 
graphs  published  in  Dow  Glycols  (2). 

Because  of  the  hygroscopic  property  of  glycols,  it  was  imperative  to 
determine  whether  changes  in  the  water  content  of  the  fluid  media  had 
occurred  between  the  times  of  utilization  of  the  fluid  media.  As  an  ex- 
ample, an  evaluation  of  the  viscosity  of  the  fluid  media  required  the 
determination  of  a  density  value  by  actual  measurement.  This,  in  turn, 
provided  a  step  toward  the  determination  of  the  percentage  of  pure  tri- 
ethylene glycol  in  the  fluid  media.  The  density  of  the  fluid  media  was 
obtained  by  the  utilization  of  a  volumetric  flask,  scales,  and  a  properly 
selected  sample.  From  this  density  and  the  density  of  pure  water  at  the 
same  temperature,  25°C,  the  specific  gravity  of  the  fluid  media  at  25°C 
was  calculated,  With  this  information  the  percentage  of  Triethylene 
Glycol  was  found  from  a  plot  of  specific  gravity,  25°C/25°Cj versus  per 
cent  glycol  by  weight  (2).  In  turn  the  viscosity  at  any  desired  temperature 
was  found  from  a  plot  of  absolute  viscosities  of  aqeous  Triethylene  Glycol 
solutions  versus  Temperature  (2).  It  was  found  that  over  a  twelve  month 
period  with  the  proper  protection  of  the  fluid  from  the  atmosphere  no 
appreciable  change  occurred  in  the  water  content  of  the  approximately 
200  lbs.  of  fluid.  It  remained  approximately  100$  pure  Triethylene  Glycol. 
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TABLE  OP  NOMENCLATURE 

o 

1.  s      =  mean  square  deviation  from  regression,, 

y°X  — 

2.  s      =  sample  standard  deviation  from  regression  =  is    . 

3.  s2   Q  -  s2    (P  ). 

y.x6     y.x  v  e 

4.  8  =  number  of  motion  picture  frames  through  which  disk  is 

measured. 

5.  Wa  =   dynamic  stability  based  on  s     =  —2 — ' . 

6.  <|^s2    -  dynamic  stability  =|^s   x  . 

7.  U*artta  =  dynamic  stability  based  on  area  between  least  squares 
'       curve  and  path  of  disk  descent,  per  unit  time. 


3.    f  =  dynamic  unstability,   tj-? ,  /j^— 


AREA 


d    -  absolute  deviations  from  regression  =  sy.x  (6-2) 
y.x 


11.  tfl  =  time  required  for  disk  to  descend  9  frames. 

12.  L_  =  measured  height  of  top  reference  line9  relative  to  base 

line. 

lj.  Lp  m   measured  height  of  disk,  relative  to  base  line. 

14.  L   =  measured  height  of  bottom  reference  line,  relative  to 

base  line. 

15.  T.L.  m        100    ■■■-   factor  for  adjusting  relative  values  to 

L„  -  L.    true  values. 

16.  h   =  (L  -  L  )  (    100   )  m   true  height  of  particle  above 

1  L   L_  -  L.      reference  line. 

17.  n  m   number  of  sample  under  statistical  consideration. 

18.  k  ■  y  intercept  of  least  squares  line. 

19.  j  =  slope  of  least  squares  line. 

20.  I  =  moment  of  inertia  of  disk. 
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TABLE  OF  NOMENCLATURE  (oonol.) 

21.  Re  =  Reynolds  number  "based  on  J  »  -jp-  . 

22.  D  =  diameter  of  disk  in  feet. 

2 J.  Y  =  average  velocity  of  disk  in  feet  per  second. 

24.  p    =  density  of  disk  in  pounds/feet  cube. 

25.  M    m     viscosity  in  British  viscosity  units. 

D  V© 

26.  Re  =  Reynolds  number  based  on  I)  ■  s  r  . 


H 

27.  D   ■  diameter  of  sphere  having  same  volume  as  disk. 

28.  t„  __   =  student's  t  for  95  per  cent  confidence  interval  =  2.571. 

U.U5 

29.  5  =  conversion  factor  of  frames  to  seconds. 

30.  x     =  time  axis  ""j 

>  on  plot  of  disk  path  (PLATE  XI?) . 

31.  y  -   height  axisj 
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TABLE  OF  DEFINITIONS 

1.  Drop  referred  to  a  sum  of  single  descents  for  all  disks  Investigated. 

2.  Descent  referred  to  a  movement  downward  of  a  single  disk. 

3.  Dynamic  stability  referred  to  the  stability  of  a  disk  during  free 
fall  in  a  container  of  infinite  proportions. 

U.       Parallax  was  defined  as  a  distortion  of  the  measured  values  of  the 
differences  between  consecutive  positions  of  a  descending  disk. 

5.  Consecutive  positions  referred  to  consecutive  frames  on  the  film  of 
the  descent  of  a  disk  where  these  positions,  or  frames,  were  separated 
by  a  constant  time  increment. 

6.  Confidence  intervals  were  defined  as  the  limits  between  which  a  true 
value  had  a  certain  probability  of  existing. 

7.  Dark  room  was  the  room  containing  the  equipment  for  measuring  disk 
heights  on  the  films  of  the  descent. 

8.  Least  squares  line  was  the  best  possible  straight  line  through  a  plot 
of  the  path  of  the  particle  descent  which  minimizes  the  squares  of  the 
deviations. 

9.  Measuring  error  referred  to  the  error  resulting  from  the  measurements 
of  disk  height  made  on  the  films  of  the  disk  descent. 

10.  Extraneous  disturbances  referred  to  disturbances  caused  by  movements 
outside  of  the  measuring  area. 

11.  Stable  Motion  (See  PLATE  XVII). 

12.  Wave  I   (See  PLATE  XVII). 

13.  Wave  II  (See  PLATE  XVIII)' 
U.  Wave  III  (See  PLATE  XVIII). 


Ill 


TABLE  OF  DEFINITIONS  (oonol.) 

15.  Oscillation  (See  PLATE  XVIV). 

16.  Spiral  (See  PLATE  XVIV). 

17.  Agile  plastic  stars  were  star  shaped  pieces  of  plastic  approximately 
1-J-  inches  in  diameter  and  -£•  inch  thick  used  normally  to  minimize 
evaporation  from  a  liquid  surface. 

18.  Moment  of  Inertia  referred  to  integral  from  center  of  the  core  to 
outside  of  the  core  of  the  differential  mass  times  the  arm  squared 
plus  integral  from  outside  of  core  to  outside  of  rim  of  the  differ- 
ential mass  times  arm  squared,  where  the  arm  is  the  distance  from  the 
center  of  the  disk  to  each  differential  mass. 
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SOAP 

Output 

Of  A 

FORTRANSIT 

Program 

For 


IPs  And  Least  Squares  Line 


Programmed 
by 

Mr.  Thomas  L.  Hamilton 

IBM  650 
Kansas  State  University 


Program  is  preceded  by  a  SOAP  package 
deck  from  FORTRANSIT  II. 
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EXPLANATION  OF  PLATE  XXVI 

Photograph  of  drop  column,  corresponding 
to  drawing  of  drop  column  presented  in 
Plate  II. 
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ABSTRACT 

This  study  pertained  to  the  characteristics  of  the  motions  of  particles 
falling  freely  in  a  fluid. 

There  were  four  objectives  to  this  investigation.  The  first  was  to 
refine  the  mathematical  definition  of  a  property  of  motion  of  each  particle 
as  it  freely  descended  in  an  infinite  expanse  of  fluid.  This  property  was 
to  be  named  dynamic  stability.  The  second  objective  was  to  actually  measure 
the  behavior  of  the  motion  of  each  particle  and  subsequently  calculate 
values  of  dynamic  stability.  The  third  was  to  observe  and  define  the  type 
of  motion  exhibited  by  each  particle  falling  freely  in  a  fluid.  The  fourth 
was  to  correlate  the  results  of  two  and  three. 

The  dynamic  stabilities  were  found  by  dropping  a  series  of  particles, 
one  at  a  time,  through  a  glass  column  containing  triethylene  glycol.  The 
descending  particles  were  filmed  on  moving  pictures  and,  from  the  individ- 
ual frames  of  this  film,  measurements  were  made  of  the  heights  of  the 
particles  at  equal  increments  of  time.  From  these  measurements  a  least 
squares  line  was  calculated  for  the  path  of  each  particle.  The  mean 
deviation  of  the  actual  path  of  the  particle  from  the  least  squares  line 
gave  dynamic  unstability;  the  reciprocal  of  this  quantity  was  taken  as 
dynamic  stability. 

These  investigations  were  conducted  with  disks  -  approximately  one 
inch  in  diameter,  a  quarter  inch  thick,  and  of  various  weights  and 
moments  of  inertia  -  descending  in  triethylene  glycol  at  88  to  90T"  at 
Reynolds  number,  based  on  Dg,  ranging  from  78  to  64.0. 

It  was  found  that  a  disk,  in  a  free-fall  descent  through  triethylene 
glycol,  exhibited  one  of  three  postulated  types  of  motion;  i.  e.  stable, 


wave  I,  and  wave  III.  It  was  further  postulated  that  wave  I  motion 
consisted  of  subtypes;  wave  I-A  and  wave  I-B.  The  distinction  between 
these  subtypes  was  shown  by  a  fluctuation  in  the  values  of  dynamic 
stability. 

A  relationship  was  found  between  all  the  postulated  types  of 
motion  in  triethylene  glycol  and  the  dynamic  stabilities, OTs  _  . 
Disk  6,  exhibiting  stable  motion,  had  dynamic  stabilities  between 
11.4-9  to  14-. 92  reciprocal  cm.  Disks  postulated  as  exhibiting  wave 
I-A,  I-B,  and  III  motions  had  ranges  of  dynamic  stabilities  of  6.02 
to  8.55  reciprocal  em.,  4.15  to  5.46  reciprocal  cm.,  and  0.44  to  2.01 
reciprocal  cm.,  respectively.  Each  range  of  dynamic  stabilities,  corres- 
ponding to  a  particular  type  of  subtype  of  motion,  was  distinct  and 
separate. 

From  a  comparison  of  the  physical  properties  of  the  disks  with 
the  types  of  motion  assigned  each  disk  it  was  found  that  a  low  dynamic 
stability  was  associated  with  a  high  moment  of  inertia  relative  to 
weight  of  the  disk.  From  a  comparison  of  Reynolds  numbers  for  each 
disk  with  the  types  of  motion  it  was  found  that  a  Reynolds  number, 
based  on  Ds,  of  less  than  one  hundred  was  associated  with  stable  motion. 

This  study  shows  that  dynamic  stability  of  disks  is  measurable  and 
therefore  must  exist.  The  measurable  dynamic  stabilities  also  may  be 
utilized  to  predict  the  types  of  motion  of  disks  descending  in  a  fluid 
media. 


